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ABSTRACT 

A catalog of Cepheid variables is used to probe the kinematics of the Galactic disk. 
Radial velocities are measured for eight distant Cepheids toward i = 300°; these new 
Cepheids provide a particularly good constraint on the distance to the Galactic center, 
Rq. We model the disk with both an axisymmetric rotation curve and one with a weak 
elliptical component, and find evidence for an ellipticity of 0.043 db 0.016 near the Sun. 
Using these models, we derive Ro = 7.66 ± 0.32 kpc and v c i rc = 237 ± 12 km s _1 . 
The distance to the Galactic center agrees well with recent determinations from the 
distribution of RR Lyrae variables, and disfavors most models with large ellipticities 
at the solar orbit. 

Subject headings: Cepheids — Galaxy: fundamental parameters — Galaxy: kinematics 
and dynamics — Galaxy: structure — distance scale — techniques: radial velocities 

1. Introduction 

The first use of Cepheids to measure kinematic parameters of the rotation curve was by Joy 
(1939), who found a distance to the Galactic center of 10 kpc; this distance was inferred from the 
measured shape of the rotation curve assuming a simple model. Cepheids have several significant 
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advantages over other stellar tracers for determining large-scale Galactic kinematics. They are 
intrinsically bright {My — —4.1), which in conjunction with variability makes them relatively 
easy to locate at large distances. Perhaps the most significant advantage of Cepheids is that 
distances to them can be determined extremely well: modern calibrations of the period-luminosity 
relation in the near infrared yield uncertainties of < 5% (Madore & Freedman 1991). Distances 



to Cepheids in the Galactic disk can be best obtained from near-infrared photometry, due both to 
the smaller PL relation scatter and because of heavy extinction by dust in the Galactic plane. 

Further studies were carried out by Stibbs (1956), Kraft & Schmidt (1963), and Feast (1967), 
using additional data on Cepheids and incorporating the use of individual reddenings in distance 
measurements. Their analysis was ultimately limited, however, by the small amount of available 
data, particularly on distant, faint Cepheids. Caldwell & Coulson (1987, hereafter CC) made an 
extensive compilation of available Cepheid photometry and radial velocities, and used the data 
in an axisymmetric rotation curve model to determine, among other parameters, the distance to 
the Galactic center (Rq). Though a significant improvement over earlier work, their models were 
also limited by the available data: CC lamented that many distant Cepheids lacked good radial 
velocities, and there were few Cepheids known at large distances from the sun, particularly toward 
directions which provide the best constraints on Rq. 



More recently, new radial velocities for many distant Cepheids have been measured ( Metzger 



et al. 19911 [Metzger, Caldwell, fc Schechter 1992] , Pont et al. 1994b), and used in models by 



Caldwell et al. (1992) and Pont, Mayor, & Burki (1994a, hereafter PMB). Uncertainties in the 
parameters of the axisymmetric rotation curve models employed, including Rq, were significantly 
improved. Under the assumptions made by the models, it is one of the most precise ways to 



measure Rq ( |Reid 1993| ) 



Many studies suggest, however, that the simple axisymmetric picture of Galactic rotation 
may not be correct. Much evidence leads to the conclusion that the Milky Way is a barred spiral 
(e.g. [Blitz fc Spergel 1991b| , [Binney et al. 1991] , [Weinberg 1992j , |Dwek et al. 1995[ ), and there are 



suggestions that even orbits near the Sun are significantly elliptical (e.g. piitz &: Spcrgcl 1991a , 



[Kuijken Sz Tremaine 1994| ). A curious result of many rotation models using Cepheids, first noted 
by Joy (1939), has been an apparent constant velocity offset of ~ 3 km s _1 between mean Cepheid 
velocities and the local standard of rest (LSR). Though the Cepheid sample has grown and 
measurements have improved for recent studies, this term has persisted. Some explanations for 
this effect have included a possible error in radial velocity zero-point or an artifact of measuring 
velocities in a pulsating atmosphere (CC). One of the results of the PMB study was to suggest 
that instead this velocity offset might be a kinematic effect explained by a bar-driven ellipticity of 
the solar orbit. 

In this paper, we present new data on eight distant Milky Way Cepheids discovered by 
Caldwell, Keane, & Schechter (1991). These Cepheids are located near the solar circle, which gives 
them high leverage in determining Rq and helps to decouple the measurement of Rq from the local 
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rotation speed. We combine the new data with a catalog of Galactic Cepheid data, and measure 
new Galactic rotation parameters using an axisymmetric model. We also study the kinematics 
using the non-axisymmetric models of Kuijken & Tremaine (1994), and measure a component of 
the local potential ellipticity. 



2. Velocities for New Cepheids 

Caldwell, Keane, & Schechter (1991) conducted a search for distant Milky Way Cepheids in 
an area near t = 300°, 6 = 0° covering 9.4 square degrees. Prom over 2000 variable stars identified 
in the survey, 37 were chosen as promising Cepheid candidates. The candidates were selected 
based on sparse I-band data, which could not provide a firm classification of the candidates. 
Additional multi-band photometry of the candidates was obtained to help confirm the identity of 



these stars ( Avruch 1991 ) 



We obtained spectra for many of the candidates in February 1991 at Las Campanas, both to 
help establish these candidates as Cepheids and to obtain radial velocities for use in kinematic 
models. To select the most promising candidates for frequent observation, we made use of V- 
and I-band followup photometry of Avruch (1991). I-band data from the original survey were 
combined with new observations to generate a detailed light curve, and new V-band observations 
were used to provide color information at different pulsational phases. Since Cepheids have a 
characteristic (V—I) color change over the course of pulsation (e.g. [Moffett fc Barnes 1985| ), we 
ranked the candidates for observation based on the slope of dV/dl as well as the appearance of 
the I-band light curve. 



2.1. Observations and Data Reduction 

Spectra of the candidates were taken with the Modular Spectrograph on the DuPont 2.5m 
telescope at Las Campanas on the nights of 25 February through 2 March 1991. We used the 
spectrograph in a cross-dispersed mode with a 150 £/mm immersion grating and a 300 d./mia 
grism cross-disperser, projected onto a TI 800x800 CCD. The primary disperser was adjusted to 
place orders 14-25 onto the CCD, providing coverage from 5000-8700 A. A 1.0 arcsec slit was used 
throughout the run, which projected to 2.2 pixels on the detector at 8400 A and gave an effective 
resolution of 60 km s _1 . Calibration frames were taken after each stellar spectrum using He-Ne 
and Fe-Ar lamps. 

The data reduction was conducted using a slightly modified version of the procedure described 
by Metzger et al. (1991). Each spectrum was flattened using an incandescent lamp exposure, and 
strong cosmic ray events were removed. Calibration lines from the associated lamp exposures were 
identified and centroided using a modified version of DoPHOT ( Mateo &: Schechter 1989] ), and fit 
across orders with a fifth-order 2-dimensional Legendre polynomial. The high-order coefficients 
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were fixed using a long lamp exposure that yielded over 300 identified lines, and the 4 lowest-order 
coefficients were fit to each calibration frame, which typically had 100 available lines. Each stellar 
spectrum was rebinned in log-A according to the calibration, and each order was separately 
extracted and sky-subtracted. 

During the course of the reduction, we noticed that the spectrum shifted abruptly between 
two positions on the chip throughout the observing run. The shift was aligned in the direction of 
the cross-dispersion, and it is possible that the grism was not well secured, flopping between two 
positions. Figure |l] shows the position of the spectrum on the chip as a function of telescope hour 
angle, from which it is clear that the flop occurs near the meridian. To account for the shift in the 
data reduction, each spectrum and its associated calibration spectrum was classified into a "high" 
or "low" group. Separate flats and high-order wavelength calibrations were made for each group 
and the associated frames were reduced within its group. Three frames that suffered a shift during 
an exposure were discarded. 

A velocity for each spectrum was calculated relative to a high signal-to-noise spectrum of 
the star HD 83443 using the Fourier quotient technique of Sargent et al. (1977). Several of the 
Cepheid candidates were too heavily reddened to provide an adequate signal for radial velocities 
using the blue orders, so we decided to use the Ca triplet (8498, 8542, 8662 A) in order 14 to 



measure individual velocities. Several CORAVEL faint southern radial velocity standards ( Maurice 



et al. 1984| ) were observed throughout the run, and were used to calibrate the effective velocity 
of the template spectrum. The individual measurements of the standards (given as an effective 
template velocity, computed using radial velocity of the standard star, as by |Metzger et al. 199l| ) 
are shown in Figure |2[ The open and filled points correspond to "low" and "high" position spectra, 
respectively. The mean template velocities of the two groups are not significantly different, 
confirming that the two groups of spectra have been successfully referenced to the same zero 
point. Table p] shows the mean measured velocities for each radial velocity standard, along with 
velocities for two metal-weak subdwarfs observed during the run. Each velocity has been adjusted 
upward by 0.4 km s" 1 to bring these velocities to a minor planet-based zero point (pVlayor 1985] ; 



Metzger, Caldwell, fc Schechter 1992] , hereafter MCS). Two error estimates are given in Table % 



£f, the formal error in the velocity from the Fourier quotient, and ey = a v /\/n, the standard error 
of the mean of the individual measurements. 

The individual radial velocities for the stars confirmed as Cepheids are given in Table H The 
other variables did not have spectra consistent with that of a Cepheid; most of these did not have 
significant Ca triplet absorption lines and thus did not yield radial velocities. 



2.2. Gamma Velocities 



Gamma velocities were computed for each of the Cepheids according to the method described 
by MCS. For each star, the radial velocity points are folded at the pulsation period and fit to a 
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"typical" radial velocity curve. The shape of this velocity curve is fixed by the photometrically 
determined period of the Cepheid, and is generated using low-order Fourier coefficients that are 
functions of the period. The periods used for these stars were taken from Avruch (1991), and 
improved with additional data below. The gamma velocity and phase are then fit using a \ 2 
minimization procedure. 

The radial velocities for each star and the curves fit to the points are shown in Figure ||. 
Fit gamma velocities for the eight Cepheids are given in Table ||, along with the formal error of 
the fit and the reduced x 2 - An additional error estimate, based on a Monte Carlo simulation, 
is also given in Table ||. This estimate takes the variation in the shape of the fit curve into 
account as well as velocity measurement error. One star thought to be a Cepheid from its light 
curve (11582-6204), but suspect due to its near-infrared photometry ( |Schechter et al. 1992 ), has a 



spectrum inconsistent with that of a Cepheid and was discarded from the candidate list. 



2.3. Infrared Magnitudes 

Mean (K) magnitudes of the new Cepheids were computed using the data of Schechter et al. 
(1992). The values of K they report are a straight average of their individual measurements; if 
the measured points are not well-spaced in phase, such an average can be biased with respect to 
the true (K). Though we expect this difference to be small given the pulsation amplitude at K, 
to obtain a slightly more accurate average we fitted a sine function to the K points using periods 
computed above. The amplitude of the sine function was scaled from the V light curve amplitude, 
computed from the data of Avruch (1991) and LeDell (1993), using the relation of Welch et al. 
(1984): Amp(K) = (0.30 ± 0.03) x Amp(V). The results of this procedure, along with formal 
errors assuming gk = 0.02 mag for each observation, are given in Table ||. Adding higher-order 
terms to the light curves from the Fourier decompositions of Laney & Stobie (1993a) had no 
significant effect on the computed (K). The star 13240-6245 has a high covariance (r ~ 0.6) 
between (K) and the epoch, largely a result of poor phase coverage. 

Also shown in Table || are the amplitudes (peak-to-peak) of the F-band light curves, epochs 
of maximum light in both V and K, and improved period estimates. The isT-band epochs were 
determined from the fit light curve, and the improved periods are selected from values listed 
by Avruch (1991) (he gives several due to the possibility of aliasing) that are most consistent 
with the K data. We find that the V maximum light lags that in K by ~ 0.27 cycles, in rough 
agreement with Welch et al. (1984). We also note that the star 13323-6224 is peculiar in that it 
has a significantly smaller amplitude than expected given its period. Overall, we obtain a tight 
formal error on the (K) magnitudes; in particular, the uncertainties are smaller than the scatter in 
the PL-K relation and hence sufficient for our purposes. (It is interesting to note that only two of 
the average magnitudes are significantly different from the straight means computed by Schechter 
et al. [1992].) It was not necessary to phase and fit the {H — K) colors, since they do not change 
appreciably over the pulsation cycle. 
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We used the E(H — K) color excesses of Schechter et al. (1992) to compute the extinction in 
K, Ak = K — Kq, for the newly discovered Cepheids. These were derived assuming an intrinsic 
color locus in the H-K/P plane: (H - K) = 0.068 + 0.024(logP - 1) (their equation 1). We 
adopted the same extinction law used by Schechter et al. for the total-to-selective extinction, that 
given by Cohen et al. (1981): Ak = \.2>9E{H — K). This can be compared with coefficients found 
in other sources: 1.7 (McGonegal et al. 1983, CIT system), 1.5 ([Clayton, Cardelli, fe Mathis 1989j 
Johnson system), 1.8 ( Rleke Lebofsky 1985| ), and 1.6 (Laney & Stobie 1993a, Carter system). 
The Caldwell et al. Cepheids have an average E(H — K) of 0.27; if we were to simply replace our 
reddening law with the average of the Ak/E{H — K) values listed (= 1.6), the result would be an 
increase in mean distance modulus by 0.06 mag for the stars in the sample. This would not be 
strictly correct, however, as the values are based on magnitudes of different systems; we use the 
Cohen et al. (1981) value, keeping in mind a possible systematic offset. 



3. Galactic Cepheid Data 

The new Cepheid data from §^ were combined with a catalog of Cepheid photometry and 
radial velocities compiled from several sources. We started with the compilation of 184 stars by 
CC. Reddenings for many additional Cepheids were obtained from Fernie (1990, 1994), and new 
radial velocities from Moffett & Barnes (1987), Mermilliod et al. (1987), and MCS. Finally, mean 
B and V photometry and many new radial velocities were obtained from the compilation of PMB. 
Known or suspected W Virginis stars (Pop. II Cepheids), based on the list of Harris (1985), were 
excluded from the sample. 

Because of the small scale height of Cepheids in the Galactic disk, the distant stars lie within 
a few degrees of the Galactic plane in projection and are thus significantly obscured by dust. To 
measure distances, accurate color excesses for the stars must be known to allow an estimate of 
the total extinction. Figure ||a shows a comparison between the color excess given by CC and the 
values of Fernie (1990). A clear trend is evident, such that the redder stars tend to have higher 
values of E(B-V) on Fernie's scale than that of CC. A linear fit to the data gives 

E B ^ V {CG) - # B _y(F90) = 0.032(±0.006) - 0.107(±0.013) x E B _ V (¥99) , (1) 

with a 0.05 mag scatter about the fit. This trend implies a difference in R = Ay jE{B — V) 
between the two systems of 10%, or about 0.3 mag in the corrected V magnitude for the more 
heavily extinguished stars (Eb — V — 1). The origin of the discrepancy is not clear, but it is 
significant: a 10% change in R corresponds to a 5% change in Rq derived from rotation curve 
models (see §|| below). A comparison between reddenings of Fernie (1990) and Dean, Warren, & 
Cousins (1978) for 94 stars in common, yields a similar but shallower slope; Fernie shows a similar 
plot for cluster stars in his Figure 1. Another comparison can be made to the recent study of 
Cepheid reddenings by Laney & Stobie (1994), who combined optical and infrared data to test the 
value of R. Figure |||b shows a comparison of their reddenings to Fernie's, revealing the same trend 
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as Figure |||a. We have adopted the same value of R as Laney & Stobie found consistent with their 
color excesses, thus to produce extinction corrections we adjust Fernie's (1990) reddenings using 
equation ID. 

The combined data set contains 294 stars (not including the eight new stars from §Q) and is 
shown in Table ||. Values reported for E(B-V) are primarily from Fernie (1990), transformed to 
our system. 

3.1. Cepheid Distance Calibration 

Distances to the Cepheids were computed primarily via the period-luminosity relation in the 
V-band (PL-V), which can be parameterized 

5 log 10 d = m — mo — a(log P — 1) . (2) 

Here m,Q corresponds to the unreddened apparent magnitude of a Cepheid with a 10-day 
period at a distance of Ro, and a is the slope of the adopted period-luminosity relation. The 
measured period and unreddened apparent magnitude for a particular Cepheid are given by P 
and m, respectively. The stars were dereddened as in CC, using the prescription for R\r derived 
from Olson (1975) and Turner (1976): 

R v = A V /E{B -V) = 3.07 + 0.28(5 - V) + 0ME(B - V) , 

whereby 

m = (V) -R v x E(B - V) . 

This extinction correction is also used by Laney & Stobie (1993), and has a slightly higher value for 
the effective Ry than the reddening laws of Savage & Mathis (1979) and Rieke & Lebofsky (1985) 
(though the last two do not give an explicit dependence on intrinsic color). The dependence of 
Ry on the intrinsic color is caused by a shift in effective wavelength of the filters when measuring 
stars of different spectral class (Olson 1975). Extinction corrections for infrared photometry are 
handled separately and are described below. 

We chose to avoid period-luminosity-color (PLC) relations and terms accounting for metallicity 
in this study for two reasons. First, the PL-V relation is thought to be only weakly sensitive to 
metallicity, both in zero-point and slope (Iben & Renzini 1984, Freedman & Madore 1990; but see 
also Caldwell & Coulson 1986), and less sensitive to metallicity than the PLC relation (Stothers 
1988). Second, the intrinsic color used in the the PLC relation is susceptible to significant error 
when this color is found by dereddening using color excesses. By making a correction based on 
color, one must make some implicit assumption about the intrinsic color of the star; any difference 
between this assumption and the actual color (perhaps due to metallicity) is amplified in the 
derived unreddened magnitude. This process can significantly increase the sensitivity of PLC 
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distances to metallicity. The uncertainties in the reddening corrections themselves tend to be 
larger than the effect of metallicity on the PL-V relation over a wide range of metal abundance 
(Stothers 1988). Further, the reddenings derived for Cepheids themselves are sometimes derived 



under the assumption of a unique color locus in the (B-V) / (V-I) plane (e.g. Dean, Warren, fc 



Cousins 1978| ). Given these sources of error, and the uncertainty in the slope of the PLC color 



term (Fernie & McGonegal 1983, Caldwell & Coulson 1986), for this study we chose to use only PL 
relations (PL-V, plus the PL-K relation for some of the models below) with no explicit correction 
for a radial metallicity gradient. The scatter in the individual distances may be slightly higher, 
but the systematic errors are easier to quantify. 

The two parameters mo and a in equation [2| determine the distances to each Cepheid in terms 
of Rq. The PL-V slope parameter a has been measured using both Magellanic Cloud and Galactic 
cluster Cepheids (Fernie & McGonegal 1983, Caldwell & Coulson 1986, CC, Madore & Freedman 
1991, Laney & Stobie 1994). Most tend to agree to within the quoted errors, and lie in the range 
—2.9 to —2.8 (with the exception of CC at —3.1). There does appear to be a significant difference 
in computed slope of the PL relation, however, depending on the period range of Cepheids used 
in the fit. While studies using open clusters to calibrate the PL relation contain data over a 
wide range of period, many exclude the longest period Cepheids from the fit as they tend to 
be somewhat brighter than an extrapolation of the PL relation of short-period Cepheids would 
indicate (e.g. Fernie & McGonegal 1983). Freedman et al. 1993 derive a separate calibration of 
the PL relation based only on Cepheids with 1.0 < logP < 1.8 to match most closely the range of 
periods in the M81 Cepheids. They find a PL-V slope of —3.35 ± 0.22, significantly steeper than 
when short-period Cepheids are used. 

A possible explanation for the discrepancy is the suggestion by Bohm-Vitense (1994) that 
most Cepheids with periods shorter than 9 days are overtone pulsators. If the short and long 
period Cepheids form two offset, steeper PL relations, then a slope measured from combining the 
two would be shallower than that measured from either set independently. More work needs to 
be done to help verify the existence of the separate PL relations, particularly in the near- infrared 
where the intrinsic scatter about the PL relation is smaller. A quick examination of the PL-K 
data of Laney & Stobie (1994) shows little evidence for short-period overtone pulsators, while 
not necessarily ruling them out. Gieren, Barnes, & Moffett (1989) find evidence against this 
hypothesis based on the continuity of BW radii across a wide range of periods. 

There is direct evidence, however, that at least some of the Cepheids in our sample are 
overtone pulsators. For example, the recently measured radial velocities for QZ Nor of —38.6 ± 0.7 
km s" 1 (MCS) and V340 Nor of -40.0 ± 0.1 km s" 1 QMermilliod et al. 1987|) confirm both as 
members of NGC 6067, and support the conclusion of Coulson & Caldwell (1985) and Moffett & 
Barnes (1986) that QZ Nor is an overtone pulsator while V340 Nor is not. Given this conclusion, 
the period-luminosity relation gives the same distance to both. Other examples of such direct 
evidence include SU Cas ( [Evans 1991 ). 
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Even with some contamination from overtone pulsators, so long as the range of periods of the 
calibrators is similar to the overall population, the derived slope and zero-point will still provide 
accurate distances (though perhaps with larger scatter). Considering that our Cepheid sample 
has a median logP ~ 0.9, we can comfortably use the shallower slopes derived from Cepheids of 
similar period, and adopt a commensurate zero point. New data on Magellanic Cloud Cepheids 
discovered in the MACHO survey ( [Alves et al. 1995| , Alcock et al. 199~5| ) will help to answer this 
question conclusively. 

The zero point of the Cepheid PL relation puts mo, the unreddened apparent magnitude 
of a Cepheid at a distance Rq, on an absolute distance scale. Various studies have yielded 
different Cepheid PL zero points, primarily due to differences in assumed extinction, metallicity or 
correction for metallicity, and the sample of stars used. Currently the most accurate methods for 
Galactic PL calibrations are those using Cepheids in clusters and associations ( [Turner 1985 ; Fernie 



fc McGonegal 1983| ), and those using the visual surface brightness (Baade-Wesselink) method 



(Gicren, Barnes, & Moffett 1989). The cluster calibrations are based on fitting main sequences for 



clusters containing Cepheids to either the Hyades or Pleiades, and the surface brightness method 
attempts to measure the radius of a Cepheid based upon accurate photometry and radial velocity 
measurements. A convenient comparison of Cepheid calibrations can be made by applying the 
calibrations to LMC Cepheids, and comparing the derived LMC distance moduli. The SMC is 
somewhat less suited to this purpose as it is thought to be significantly extended along the line 
of sight. Feast & Walker (1987) give a comprehensive review of Cepheid calibrations up to that 
time, and conclude that for a Pleiades modulus of 5.57, the LMC lies at a true distance modulus 
of 18.47 ± 0.15. This estimate is based on the same extinction law used here. More recently, using 
updated V-band data, CC determine an LMC modulus of 18.45, and Laney & Stobie (1994) find 
18.50 ± 0.07, both assuming the same Pleiades modulus and extinction law. The visual surface 
brightness calibrations currently yield Cepheid distance moduli larger by ~ 0.15 mag on average 
( Gieren &: Fouque 1993| ), and give a distance modulus for the LMC of 18.71 ± 0.10 mag. While 



significantly different, the Gieren &; Fouque data appear to have an asymmetric distribution that 
may make the distance moduli too large: the four calibrators they discard as being significantly 
discrepant all have distance moduli too large by > 0.6 mag. Though this does not completely 
resolve the discrepancy, for this study we have chosen to adopt the cluster calibrations. On this 
scale, our adopted V-band calibration is 

My = -4.10-2.87(logP- 1). (3) 

The internal uncertainty in the zero point (exclusive of any systematic error in the Pleiades/LMC 
distance) is estimated to be ~ 0.07 mag. 

To incorporate the near-infrared data on the newly-discovered Cepheids, we use the 
period-luminosity relation in the -RT-band with an appropriate calibration and extinction law. 
The calibration zero point must give distances commensurate with those derived from ^-band 
data, and thus we again normalize the zero-point to an LMC modulus of 18.50. After making 
this correction, the Ph-K calibration of Welch et al. (1987) gives Mk = —5.66 — 3.37(logP — 1), 
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with the K magnitudes on the same system (Elias et al. 1982) as the Schechter et al. (1992) 
photometry. Madore & Freedman (1991) give a self-consistent calibration based on a sample of 25 
LMC Cepheids, each with photometry in both V and K, finding Mx = —5.70 — 3.42(logP — 1), 
identical to within quoted errors. The Madore & Freedman calibration is also consistent 

with our adopted My. Laney &: Stobie (1994) give a calibration of the PL-K relation in a 
slightly different photometric system; after converting to the Elias et al. (1982) system using the 
transformation of Laney & Stobie (1993b), and adjusting to an LMC modulus of 18.50, we find 

M K = -5.70(±0.04) - 3.40(±0.05)(log P - 1) . (4) 

The scatter of the individual stars about the period-luminosity relation is significantly smaller 
in K than V, 0.16 mag rms vs. 0.25 mag rms, and hence the internal error associated with the 
zero point is correspondingly smaller at 0.04 mag. Since the quoted uncertainties in the Laney & 
Stobie (1994) calibration are the smallest of those quoted above, and that the relation is almost 
identical to the others, we adopt equation |3] for our models. 

There is a small discrepancy between the LMC moduli derived from V and K data when 
using Galactic cluster calibrations. Measurements of LMC distance modulus in the K band from 
the above references typically yield a value of 18.55-18.60, some 0.05-0.10 higher than the V 
calibration. This discrepancy could be due to a number of factors, including a difference in mean 
metallicity between Galactic and LMC Cepheids. As the bolometric PL relation for Cepheids is 
thought to be essentially independent of metallicity Qlben fc Renzini 1984 ), the zero point of the 



PL-V and PL-K relations may be expected to differ by ~ 0.03 mag over the metallicity difference 



of the Galaxy and the LMC from the variation in bolometric correction alone (Laney & Stobie 



1994| , (Stothers 19881) . Another possible source of systematic error arises from the correction for 



extinction: this is substantially larger for the Galactic calibrators, which have a mean E(B-V) of 



0.65 QFeast fc Walker 1987| ), than the LMC Cepheids, which have an E(B-V) of about 0.14 mag. 
A reasonable error of 0.1 in the adopted Ry value would thus produce an apparent distance offset 
between the two of 0.05 mag. Since the difference here is only slightly greater than 1 <r, no useful 
limits can be placed on R (or Ay — Ak). However, we discuss below some implications of the 
kinematic distance scale using the newly discovered Cepheids on the adopted reddening law. 



4. Milky Way Rotation 

Figure || shows a geometric picture of a star with angular velocity at a distance D from the 
Sun, and at a distance Rq from the center of the galaxy. The local standard of rest (LSR) rotates 
with velocity 0q, and the Sun moves with a velocity peculiar to the LSR of (uq,vo,wo) in the 
coordinates (R, 0, Z). For our models we assume that the height above the disk, \z\ = D\ sin6|, 
is sufficiently small that the potential is dominated by the disk, and therefore a thin-disk model 
adequately represents the orbits (i.e. the primarily rotational orbits of the stars are decoupled 
from their vertical motion). This condition is met by the classical Cepheid population as they 
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are confined to the disk with a scale height of 70 pc (average absolute distance from the mean 
plane jKraft fc Schmidt 1963 ). We also fix the Z-component of the Sun's motion, wq, at 7 km s _1 
( Dclhayc 1965| ) as it is neither constrained by the data nor does it affect the model. 

The rotation models we employ take the form of a mean velocity field @(R,4>) with 
components in the R and (ft directions. The mean velocity of the stars at (R,(ft), as measured from 
the Sun, is 

v* = ®r(Ro, 0) cos £cos b — @r(R, (ft) cos((ft + £) cos b 

—Q^Rq, 0) sin £ cos b + 0^>(-R, (ft) sin((ft + £) cos b — Vq* . (5) 

The Sun's peculiar motion relative to the LSR in the direction of the star, Vq*, is given by 

v ©* = — uq cos I cos b + v o sin £ cos b — wq sin b . (6) 



The data are fit to the models using a non-linear x 2 minimization program. Free parameters 
were determined by fitting the measured radial velocities for each Cepheid to model velocities 
generated from the other measured quantities (£, b, (V), P, and E(B — V)) via equation [|. Each 
measured velocity was weighted as in CC using the estimated radial velocity dispersion added in 
quadrature to the effective velocity error introduced by the distance measurement: 

af = a 2 v + al(dv r /dd) 2 . (7) 

The dispersion in the radial velocities is a combination of measurement error and the intrinsic 
velocity dispersion of the stars in the disk (the latter dominating), and was taken to be a v = 11 
km s . The random error in the distance from all sources (measurement, extinction correction, 
and PL dispersion), cr^, was assumed to be 0.2 mag. 



4.1. Axisymmetric Models 

The initial models we use to derive parameters of Galactic rotation are based on a linear, 
axisymmetric rotation curve in a manner similar to CC. An axisymmetric rotation curve is given 
by Qr = and Q^R) independent of (ft. If we make the approximation that the rotation curve is 
linear, i.e. that 

e *<*> M e » + (H " *> (§) «,= e ° + (r " 11 it] rrf 

equation || reduces to (cf. [Mihalas &: Binney 1981| , CC) 

v* = -2AR (l - sin ^ cos 6 - t; * . (8) 
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Here we have defined r = R/Rq, and A is Oort's constant 



2 



@o_(d&\ 
R \dRJ 



Ro 



The transformation of the heliocentric distance d = D/Rq to the Galactocentric distance r is given 
by 

r 2 = l + d 2 - 2dcos£. (9) 

While the expansion of the rotation curve is to first order in r, and thus is valid primarily for stars 
near the solar circle, no approximation is made for small d. 

An additional parameter was added to compensate for a possible zero-point offset in the 
radial velocities (5v r ). The full model for the measured heliocentric radial velocity of a star is then 

v r = — 2ARq ^1 ^ sin £ cos b 

+ no cos I cos b 

— vo sin i cos b 

— wo sin b — 5v r . (10) 

Note that the sign of no follows a Galactic radial convention and is opposite that used by CC. 
The model has parameters 2ARo, mo, a, no, vo, wo, and 6v r ; a was fixed according to the PL 
relation adopted, wo was fixed as described above, and the remaining parameters were fit using 
the x 2 minimization routine. 

Table § gives a summary of the results of the axisymmetric models. Model Al includes all 
of the sample Cepheids, including the 8 new Cepheids with distances determined from K-band 
photometry and reddenings. Distances to the remaining stars were computed using the color 
excess scale given by equation |l]. Model Al.l is identical except for the exclusion of 6 stars with 
large residuals in model Al: FF Car, FM Car, BB Gem, VW Pup, AA Ser, and V Vel. The model 
parameters remain essentially the same, while the associated uncertainties fall (the six excluded 
stars contribute almost 25% of the residual error in model Al). The error on Ro and 2ARq are 
independently less than 5%; the high covariance of the two errors (see below) implies an even 
tighter constraint on A = 2ARo/2Ro = 15 A ± 0.3 km s _1 kpc -1 . It is important to note that the 
errors given are internal errors (those due to scatter about the adjustable parameters), and do not 
include systematic uncertainties associated with the fixed model. Model Al.2 uses the same stars 
as Al.l, but employs the color excesses of Fernie (1990) directly, without adjustment. The most 
significant effect is a shortening of the distance scale by ~ 5%. 

As a check, we ran our modeling software on the data used by CC in their study. The 
parameters are listed as model A2 in Table and agree well with those determined by CC. The 
largest difference between the Al and A2 parameters, aside from the reduced uncertainties, is the 
value determined for 2ARq. Model A2.3 was run on the same set of Cepheids, but with updated 
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velocities from Table |5| and reddenings from Fernie (1990). Again, mo decreases by about 5%, 
while 2ARq increases by 0.8cr. Oort's A constant effectively increases from 14.2 km s _1 kpc -1 to 
15.7 km s _1 kpc -1 when the new data are used. Models A2.2 and A2.3 show the effect of adding 
the new velocities and reddenings separately. 

We also fit our models to the data set used by PMB as a further test. Model A3 uses 
their pruned set of 266 stars, and further the parameter 5v r was fixed at zero as in their study. 
Our model reproduces their results closely, though our estimates of uncertainty in the model 
parameters are somewhat larger than reported by PMB. 

Since the PMB study employed color excesses on the F90 scale, the effective distance scale 
should be the same between model A3 and A1.2. However, the computed Rq is somewhat larger 
in model A3 than Al.2, indicating a difference in the implied distance to the Galactic center. 
At first we guessed this might be due to the effect of adding the eight new distant Cepheids 
to the model-the high leverage on Rq for these stars could produce the change (particularly 
since distances were computed using a separate PL relation). Deleting these stars, however, 
only increased the uncertainty of Ro without changing the most likely value. We found that the 
discrepancy was caused by including the Sv r parameter in the model: the values determined for 
Rq and 2ARq fall by 1-cr (model A3.1 vs. model A3), and resolve the discrepancy. We conclude 
that not including the 5v r in the axisymmetric models skews PMB's Rq measurement slightly 
high. Allowance for this term can be made in non-axisymmetric models, as PMB suggest, and is 
treated in § |4.2| . 

We can further determine the effect of the new Cepheid data by examining the covariances 
between model parameters. Table |7| shows the covariances of the model parameters determined 
using different sets of data. Covariances are expressed here as correlation coefficients of the 
projected data (Bevington 1969): 

_ v 
f"ij — 

Sii Sjj 

where the sfj are elements of the covariance matrix. The reduction in the covariance between mo 
and 2ARq in the model including the new Cepheids (Al.l) over that without (A1.4) is due to the 
advantageous placement of the new Cepheids. Because they lie at a Galactocentric radius near 
the solar circle, the mean radial velocity of the stars in the model is close to zero, independent of 
rotation speed. (In an axisymmetric model, stars lying at the same radius rotate together in a 
ring with no relative velocity, save for random motion.) Thus adding only a few stars serves to 
decouple the two parameters in the model. Figure [7| shows a comparison between constant \ 2 
contours in mo, 2ARq space with and without the eight stars. 

The eight new CKS Cepheids turn out to provide a useful constraint on Rq by themselves, 
assuming axisymmetric rotation. Fitting the eight Cepheids to the model, we derive Rq = 8.1 ±0.5 
kpc, assuming vq = 14 km s —1 . Thus with just eight stars, knowing vq gives us a distance to the 
Galactic center to 6% precision. If we remove the nearest Cepheid, 13323-6224 at 0.3i?o, most of 
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the remaining covariance between Rq and 2ARq is eliminated. The error reduces to ~ 4%, with 
Rq = 7.95 ± 0.31 kpc. The parameters i?o,i>o have a high covariance in this model; decreasing 
vq by 4 km s _1 reduces the derived value of Rq by 5%. Again, all of these models rely on the 
assumption of an axisymmetric rotation curve; in the next section, we relax this assumption. 



4.2. Non- Axisymmetric Models 

The persistence of a significant Sv r in the Cepheid data, even after significant improvements 
in radial velocities and distances, leads to the conclusion that either some systematic error is 
present in measuring 7— velocities of Cepheids, or that an axisymmetric model is not sufficient to 
describe the rotation curve. As many Cepheids in open clusters now have accurate 7-velocities 



within 0.4 km s _1 of the cluster mean velocity (e.g. Mermilliod et al. 1987 ), the systematic error 
in measuring 7 is probably small. Based on a comparison with an N-body simulation, PMB 
suggest that the effect could be due to non-axisymmetric motion driven by a central bar of ~ 5 
kpc in extent. Here we examine the shape of the local rotation ellipticity directly using a simple 
non-axisymmetric model. 

Our non-axisymmetric models are based on those of Kuijken & Tremaine (1994, hereafter 
KT). The rotation curve is produced by a primarily axisymmetric potential with a small m = 2 
perturbation, with minor axis in the direction 0^. The circular velocity and potential ellipticity 
have a power-law dependence on radius: 

/ t> \ a / R \ P— 2" 

Vc{R) = @0 [T ) ; e(jR)=£o UJ ' (U) 
The rotation curve is given by equation (||) and by (KT, eqn. 5a): 

Q R (R,<P) = j3iv c (R) [s(R) cos 20 -c(R) sin 20], 

9^(i?,0) = v c (R) - f3 2 v c {R)[c{R) cos 20- s(R) sin 2(f)], (12) 



where we have defined 



and two orthogonal projections of the ellipticity 

c(R) = e(R) cos 2<j) b , s(R) = e(R) sin 2<j) b . (14) 

The c(R) and s(R) parameters correspond to components of the ellipticity that are symmetric and 
antisymmetric, respectively, about the Sun-center line = 0. 

As pointed out by KT, if the Sun lies near a symmetry axis of a non-axisymmetric distortion, 
it is not easily detected using data such as our sample where d ^ Rq. In particular, the parameters 
c and Rq in the models have a high covariance, and thus do not produce independent information. 
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If such a distortion were present (c non-zero), the derived parameters Ro and 2ARq can deviate 
significantly from their true values. To judge the size of the effect for our data, we have fit a simple 
model with a flat rotation curve and constant ellipticity (p = a = 0), and allow the symmetric 
ellipticity component c to vary along with the other parameters. For this model, the predicted 
mean radial velocities reduce to 

v r = 0o y lj sin £ cos 6 + 

+— c(2dcos t — 1) sin£ — i> * . (15) 
r 

As expected, if we fit for all six parameters, we find a covariance between 2ARq, mo, and c 
of 0.95 — too large to make any independent constraints on these parameters. We can, however, 
fix the symmetric ellipticity in the models and determine how this affects the other derived 
parameters. Figure || shows how too and 2ARq vary as a function of the ellipticity parameter. 
From a combination of non- kinematic estimates of Ro ( [Reid 1993| ), we can deduce a weak 
constraint on the ellipticity at the solar circle —0.08 < c(Rq) < 0.14. 

Because c is degenerate with other parameters, we fix it at zero and investigate the 
antisymmetric term s. For these models we incorporate the same data set used in model Al.l, 
and for the primary model assume a flat rotation curve and constant ellipticity {a = p = 0). Here 
we also assume that the velocity offset 5v r is a kinematic effect of the other model parameters, 
and therefore fix 5v r = 0. Table [8| shows the results of several models. Model Bl gives results for 
the primary elliptical model, and shows a 2.5a detection of ellipticity, s(Rq) = 0.043 ± 0.016. The 
B2 models vary the extinction coefficient R, and B3 models vary p and a; the ellipticity detection 
is robust in all except the p > models. Note that while the implied circular velocity increases 
rapidly as a increases, Oort's A remains roughly constant. Velocity residuals for model Bl are 
shown in Figure []; there is some hint that systematic velocity structure not accounted for by the 
simple models is present, perhaps due to the influence of spiral structure. 



5. Discussion 

The uncertainties quoted in Tables B48] correspond to internal error estimates in the models, 



derived from the scatter of the data about the fit model. As described in §3.1 above, systematic 
uncertainties in the distance scale are also present. To summarize, the internal uncertainty in the 
l/-band zero point is ~ 0.07 mag, to which we add an estimated 0.1 mag error in the Pleiades 
modulus; similarly, the internal uncertainty in the iC-band calibration is 0.04 mag. Thus we adopt 
a calibration error of 0.12 mag and 0.11 mag for V and K zero points, respectively. There is also 
an uncertainty associated with the reddening scale for the y-band distances; however, we note 
that since we have a good measurement of Rq based only on i^-band photometry, this uncertainty 
does not contribute to the overall distance scale error. The reddening uncertainty in the K-hscad 
corresponds to ~ 0.05 mag in Rq. The -ftT-band distance scale agrees closely with that derived 
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from the full models using the Dean et al. (1978) optical reddening scale, thus we adopt that 
scale for our derived distances. The largest systematic uncertainty in measuring Rq, however, 
derives from the inability of the rotation models to adequately constrain c(Rq). KT suggest, 
based on a combination of local kinematics and global HI distribution, that c(Rq) = 0.082 ± 0.014. 
This corresponds to overestimating Rq in the kinematic models by ~ 15% (see Figure |8|). For 
now, we keep this component of the systematic uncertainty separate; below we see that based 
on the RR Lyrae distance scale, c(Rq) is likely small. Our standard model Bl then gives 
R = 7.66 ± 0.32 ± 0.44 kpc and v c (R ) = 237 ± 12 ± 13 km s -1 . Computing Oort's A from these 
gives A = 15.5 ± 0.4 ± 1.2 km s _1 kpc -1 . 

Since the models of CC, the number of Cepheids with data adequate for use in rotation 
models has increased by 50%, in particular the addition of stars at large distance from the Sun 
that provide good leverage for measuring scale and shape parameters. The internal uncertainty 
in Rq and 2ARq has decreased by approximately a factor of two, while remaining consistent in 
their mean values. In comparison with the recent PMB models, our measurement of Rq is slightly 
smaller than that found by PMB, 8.1 kpc; much of this difference can be attributed to the absence 
of a direct ellipticity term in the PMB models. This also appears to skew the measurement of 
2ARq to the high side; our derived values differ by almost 2-cr. We find that other differences 
between our models, including a higher-order expansion of the rotation speed and a different 
minimization algorithm, produce negligible differences in the fit parameters. 

Our results compare favorably with other methods of determining Rq. Reid (1993) computes 
a "best value" for Rq, based on a weighted average of many techniques, of 8.0 ± 0.5 kpc, in 
agreement with the results from this paper. In particular, Reid divides measurements into several 
categories; our present results fall into two different categories, one that is based primarily on the 
calibration of the if-band Cepheid PL relation, and one that is based on a kinematic rotation 
model. Our distance is somewhat larger than that implied by a direct measurement of the distance 
to Sgr B2 at the Galactic center of 7.1 ± 1.5 kpc ( [Reid et al. 1988 ), though still well within the 
errors. A direct measurement of the proper motion of Sgr A* by Backer & Sramek (1987), when 
interpreted with our measured circular velocity, translates to 8.3 ±1.0 kpc, again consistent with 
our measurement, though the uncertainty is a factor of 2-3 larger. 

We can also use Rq as a standard length to compare the Cepheid and RR Lyrae distance 
scales. Based on new near-infrared photometry of RR Lyraes in Baade's window, Carney et al. 
(1995) measure Rq to a precision of about 5%. This result is much less sensitive the uncertain 
extinction correction than the previous optical work (e.g. Walker fc Terndrup 1991 ), and is also 
independent of the rotation curve shape. Carney et al. quote two different values for Rq: 7.8 ± 0.4 
kpc, based on the zero-point calibration of Carney, Storm, & Jones (1992, CSJ); and 8.9 ± 0.5 
kpc, based on an RR Lyrae calibration from the LMC corresponding to an LMC modulus of 18.50. 
As noted by Walker (1992), the Cepheid LMC distance modulus is 0.2-0.3 mag higher than that 
derived from LMC RR Lyraes on the CSJ zero point. Our result for Rq agrees quite well with the 
CSJ-scale distance, and is about 2-a below the LMC RR Lyrae-based calibration. The systematic 
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change in our Rq measurement implied by increasing c(Rq) is in the opposite sense to that implied 
by adopting an LMC RR Lyrae distance modulus of 18.50: if the KT suggestion that c(Rq) ~ 0.08 
is correct, it would only make our estimate for Rq shorter and the discrepancy with the LMC RR 
Lyrae calibration larger. Thus the CSJ calibration of RR Lyraes and the Cepheid distance scale 
appear to be commensurate to within ~ 0.15 mag, and the ellipticity term c(Ro) < 0.04. This 
implies that a small but probably real discrepancy exists between Cepheid and RR Lyrae LMC 
distances, and that the intrinsic magnitudes of the populations in the Galaxy and LMC may be 
different ( van den Bergh 1995| ; Gould 1994 ).. The two distance scales appear to be in agreement, 



however, when Cepheids are compared with globular cluster horizontal branch magnitudes in M31 



([Ajhar et al. 1996|) 



A positive value of s(Rq) such as we find indicates that the Sun's orbit is elliptical with major 
axis in the quadrant < <p < n/2, illustrated in Figure [l0| (see also Schechter 1996). PMB suggest 
that the orbital ellipticity could be due to a bar of roughly 5 kpc radius, based on results of a 
numerical simulation; this is somewhat larger than the <^ 3 kpc scale suggested by photometric 
observations (e.g. Weinberg 1992| ). The kinematic signature of the smaller-scale bar at Rq is 



predicted to be ^ 1% in the models of Weinberg (1994), thus would be a strict lower limit for the 
scale of the bar if the outward motion of the LSR is such a kinematic signature. The B3 models 
in Table || with p < correspond to these inner bar models with ellipticity falling at Rq. 

Our measurement of s(Rq) indicates that the Sun has a small outward velocity component 
with respect to the GC, of e R (R ,0) = 10.2 ± 3.8 km s" 1 for our fiducial model Bl. This model 
has p = 0, i.e. the ellipticity is constant with radius, which predicts that stars toward £ = 180° 
should have zero mean velocity. Observations of old disk stars at the anticenter, however, indicate 
a mean motion of ~ 6 km s _1 in the sense that the LSR and the outer Galaxy stars are moving 
apart (Lewis fc Freeman 1989, Mctzger fc Schechter 1994). This would suggest that s(R) increases 



with radius, such as might be the effect of a triaxial halo that gradually dominates an axisymmetric 
disk potential as R increases. Orbits at larger radius would thus have larger ellipticity, implying 
a larger < v r > (so long as a > 0), and would thus appear to have a net outward motion when 
viewed from inner orbits. This corresponds to models B3 with p > 0; when we fit to these models, 
however, the magnitude of local ellipticity shrinks and becomes negligible with p = 1. The outward 
motion is consistent with the 14 km s _1 value of the Blitz & Spergel (1991a) model derived from 
outer Galaxy gas motions; however, warps and possible lopsided structure of the outer gas make 
the velocity structure difficult to interpret ( [Kuijken 1992j ). 

One of the obstacles to drawing strong conclusions from models with different p and a is that 
the distribution of the Cepheids in the model is both asymmetric with respect to the GC and 
concentrated about R = Rq. Constraining s(R) is best accomplished when the antisymmetric 
component can be measured directly on both sides of the Galactic center. To help fill out the region 
< £ < 90°, a survey for distant Cepheids has been conducted toward I = 60° with some success 
QMetzger 199ip . Further surveys are in progress to find Cepheids at even greater distances on the 



opposite side of the galaxy; if a sample of such Cepheids were available with radial velocities, good 
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constraints could also be placed on c(R) (KT). 

It is clear that given the uncertainties in estimating total U-band extinction toward heavily 
obscured Cepheids (§|3|), the future direction of Cepheid kinematic models will be toward using 
infrared photometry. This will be particularly important for new samples of distant stars; as 
an example, some of the newly-discovered Cepheids toward I = 60° have E(B-V) of ~ 4 mag. 
An uncertainty in Ry of only 0.1 corresponds to a minimum distance error of 20%, while in K 
the corresponding extinction and error is reduced by a factor of 10. Future surveys for distant 
Cepheids, particularly those to find Cepheids at distances > Ro, are also best conducted in the 
infrared. The small scale height of Cepheids is of order that of the gas and dust in the disk, so the 
total extinction to such stars on the opposite side of the Galaxy is prohibitive in the optical. 

6. Summary 

We have obtained new 7- velocities for the eight Cepheids discovered by CKS, and combine 
these with .fT-band photometry of Schechter et al. 1992 for use in disk kinematic models. Using 
a catalog of Cepheid magnitudes, color excesses, and radial velocities collected from many 
sources, we examine Galactic rotation parameters in the context of both axisymmetric and 
non-axisymmetric models. In our adopted model, which assumes a flat rotation curve and constant 
ellipticity near the sun, we find 

R = 7.66 ± 0.32 ± 0.44 kpc, 
v c (R ) = 237 ± 12 ± 13 km s" 1 , 
s(R Q ) = 0.043 ±0.016, 
A = 15.5 ± 0.4 ± 1.2 km s" 1 kpc" 1 . 

The errors quoted are internal and systematic, respectively. An additional systematic error is also 
present from the unknown component of ellipticity that is symmetric about the sun-center line, 
but comparisons with recent measurements of Ro from RR Lyraes indicate that this component 
should be small. An estimate of Ro can also be made directly from 7 Cepheids near the solar 
orbit with infrared photometry, independent of the rotation speed, giving 7.95 ± 0.31 kpc for 
zero ellipticity and 7.61 ± 0.30 kpc for 4% ellipticity. There is good agreement between the new 
values and some previously published estimates, with significantly smaller uncertainty in the 
present results. Future progress in determining rotation curve scale and shape will rely on infrared 
searches and photometry, both to find more distant stars through heavy extinction and to measure 
the distances more accurately. 
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Table |l|. Radial Velocity Standards, Feb 91 



Star 


V r 


£ f 


e v 


n 


HD 24331 


26.6 


1.7 


1.3 


3 


HD 39194 


16.4 


1.7 


0.8 


4 


HD 48381 


40.1 


1.5 


0.5 


4 


HD 83443 


27.4 


1.2 


0.7 


6 


HD 83516 


43.3 


1.2 


1.5 


7 


HD 101266 


21.5 


1.3 


1.5 


6 


HD 111417 


-18.9 


1.2 


0.8 


7 


HD 176047 


-41.5 


1.3 


0.6 


5 


CPD -43° 2527 


19.3 


1.3 


1.2 


5 


HD 74000 


205.5 


4.6 


1.6 


4 


HD 140283 


-171.2 


3.2 


3.6 


3 
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Table Cepheid Radial Velocities, Feb 91 



JD 


V r 


<7 


JD 


V r 


a 


-2400000 


km s _1 




-2400000 


km s _1 




11447-6153 






12003-6213 






48312.748 


23.8 


4.0 


48312.730 


-3.9 


4.1 


48313.649 


31.0 


3.6 


48313.754 


-12.3 


3.8 


48314.661 


42.8 


3.2 


48314.770 


-3.5 


3.4 


48315.651 


52.0 


4.3 


48315.728 


11.8 


4.1 


48316.726 


11.3 


5.3 


48316.806 


24.5 


3.4 


48317.612 


13.9 


5.2 


48317.709 


20.0 


4.6 


11465-6209 






13190-6235 






48312.679 


-12.4 


5.7 


48312.794 


-35.4 


3.8 


48313.581 


-13.9 


4.6 


48313.812 


-22.6 


4.0 


48314.598 


-2.7 


3.0 


48314.820 


-14.2 


3.1 


48315.606 


-15.4 


4.4 


48315.817 


-8.5 


4.3 


48316.676 


-18.5 


4.0 


48316.877 


-24.5 


4.3 


48317.577 


-0.4 


3.5 


48317.770 


-45.4 


4.7 


11492-6257 






13240-6245 






48312.847 


10.2 


7.3 


48313.765 


-7.5 


3.7 


48313.682 


-6.1 


4.3 


48314.831 


-9.4 


3.8 


48314.742 


4.8 


5.6 


48315.880 


-26.6 


5.0 


48315.712 


20.0 


6.6 


48316.830 


-32.2 


5.7 


48316.757 


-6.7 


5.2 








48317.676 


-3.1 


4.6 








11521-6200 






13323-6224 






48312.693 


31.4 


4.4 


48313.782 


-39.7 


4.8 


48313.596 


4.4 


3.9 


48314.851 


-54.9 


4.9 


48314.695 


12.6 


3.6 


48315.869 


-53.3 


4.5 


48315.620 


14.1 


4.9 


48316.847 


-41.9 


3.7 


48316.701 


19.9 


3.1 


48317.799 


-37.2 


4.5 


48317.638 


26.0 


3.9 
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Table |3[ Gamma Velocities of New Cepheids 



Cepheid 


7 


Vfit 


&MC 




Period 


11447-6153 


29.1 


2.0 


1.9 


1.5 


6.4282 


11465-6209 


0.9: 


4.0 


3.8 


4.4 


11.0984 


11492-6257 


6.7 


1.6 


2.5 


0.5 


3.6798 


11521-6200 


20.4 


1.7 


1.7 


1.2 


6.6039 


12003-6213 


2.3 


1.8 


2.2 


1.1 


9.1266 


13190-6235 


-31.7 


2.2 


3.2 


1.9 


10.1576 


13240-6245 


-29.9: 


3.3 


8.6 


0.6 


15.0598 


13323-6224 


-48.3 


3.0 


2.1 


2.2 


4.2424 



a Avruch (1991). 



Table ^. Calculated Parameters for New Cepheids 



Cepheid 


Amply) 


(K) 




° Vmax 


° Kmax 


Period 


11447-6153 


0.65 


10.092 


0.013 


7962.97 


8337.43 


6.4282 


11465-6209 


0.88 


8.558 


0.014 


7968.7 


8337.84 


11.0984 


11492-6257 


0.62 


10.649 


0.009 


7959.80 


8336.00 


3.6529 


11521-6200 


0.63 


9.645 


0.014 


7958.40 


8335.16 


6.5763 


12003-6213 


0.72 


9.125 


0.012 


7958.67 


8339.86 


9.0131 


13190-6235 


0.82 


8.739 


0.015 


7958.76 


8342.12 


10.3001 


13240-6245 


1.21 


7.618 


0.022 


7966.94 


8337.60 


15.2158 


13323-6224 


0.46 


7.910 


0.012 


7961.75 


8336.63 


4.2424 



Modulo 2,440,000; a ~ 0.025 x period. 
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Table g. Cepheid Data 
Cepheid I b v r a log P (V) b > c E(B-V) 



AQL 


ETA 


40. 


.90 


-13 


,10 


-14. 


.8 


0, 


.8559 


3. 


.897 


0, 


.164 


AQL 


U 


30. 


,90 


-11, 


,60 


1. 


. 1 


0. 


,8466 


6. 


.448 


0. 


.389 


AQL 


sz 


35, 


.60 


-2. 


.30 


10. 


,5 


1, 


,2340 


8. 


,617 


0, 


.607 


AQL 


TT 


36. 


.00 


-3. 


.10 


1. 


,9 


1. 


.1384 


7. 


.147 


0. 


.476 


AQL 


FF 


49, 


,20 


6. 


.40 


-17. 


,4 


0. 


,6504 


5. 


,372 


0. 


.232 


AQL 


FM 


44, 


,30 


0. 


.90 


-7. 


,0 


0. 


,7863 


8. 


.268 


0. 


.611 


AQL 


FN 


38, 


,50 


-3. 


.10 


13. 


, 1 


0. 


.9769 


8. 


.382 


0. 


.489 


AQL 


V336 


34. 


.20 


-2. 


.10 


11. 


,5 


0. 


,8636 


9. 


.875 


0. 


.610 


AQL 


V496 


28. 


,20 


-7. 


.10 


7. 


.9 





,8330 


7. 


,720 


0. 


.402 


AQL 


V600 


43, 


,90 


-2. 


.60 


3. 


, 1 


0. 


,8597 


10. 


,034 


0. 


.812 


ARA 


V340 


335, 


,20 


-3. 


.70 


-82. 


.2 


1, 


,3183 


10. 


,242 


0. 


.547 


AUR 


Y 


166. 


,80 


4. 


.30 


8. 


,5 


0. 


.5865 


9. 


.614 


0. 


.385 


AUR 


RT 


183. 


,10 


8. 


.90 


21. 


,0 


0. 


.5715 


5. 


,450 


0. 


.076 


AUR 


RX 


165, 


,80 


-1. 


.30 


-23. 


,3 


1. 


,0654 


7. 


.664 


0. 


.278 


AUR 


SY 


164, 


,70 


2. 


.10 


-3. 


,5 


1. 


,0062 


9. 


.081 


0. 


.439 


CMA 


RW 


232, 


.00 


-3. 


.80 


50. 


.0 





,7581 


11. 


.146 


0. 


,520 


CMA 


RY 


226. 


.00 


0. 


.30 


32. 


.9 





,6701 


8. 


,105 


0. 


.253 


CMA 


RZ 


231. 


.20 


-1. 


.10 


24. 


.6 





,6289 


9. 


,698 


0. 


.448 


CMA 


ss 


239, 


,20 


-4. 


.20 


73. 


. 1 


1, 


,0921 


9. 


,939 


0. 


.524 


CMA 


TV 


227. 


,20 


-2. 


.40 


39. 


.0 


0. 


,6693 


10. 


,566 


0. 


,555 


CMA 


TW 


229. 


.20 


0. 


.10 


66. 


,5 


0. 


.8448 


9. 


.557 


0. 


.351 


CMA 


vz 


239. 


.90 


-4. 


.40 


40. 


, 1 


0. 


.4950 


9. 


,371 


0. 


.466 


CAM 


RW 


144, 


,90 


3. 


.80 


-26. 


.5 


1. 


.2152 


8. 


,657 


0. 


.614 


CAM 


RX 


145, 


,90 


4. 


.70 


-36. 


.2 


0. 


.8983 


7. 


,685 


0. 


.542 


CAR 


1 


283. 


.20 


-7. 


.00 


3. 


.6 


1 


,5507 


3. 


.735 


0. 


.183 


CAR 


U 


289. 


.10 


0. 


.10 


1. 


,7 


1. 


,5883 


6. 


,281 


0. 


.285 


CAR 


V 


275, 


,30 


-12. 


,30 


13. 


.9 


0. 


,8258 


7. 


,375 


0. 


.187 


CAR 


Y 


285, 


,70 


-0. 


.30 


-14. 


,5 


0. 


,5611 


8. 


,102 


0. 


.190 


CAR 


SX 


286, 


,70 


1. 


,30 


-11. 


,4 


0. 


,6866 


9. 


,100 


0. 


,323 


CAR 


UW 


285. 


,60 


-1. 


.80 


-14. 


,4 





.7280 


9. 


,441 


0. 


.441 


CAR 


UX 


284. 


,80 


0. 


.20 


8. 


, 1 


0. 


.5661 


8. 


,266 


0. 


.141 


CAR 


UY 


287, 


,20 


-3. 


.20 


4. 


,0 


0. 


,7438 


8. 


,928 


0. 


.200 


CAR 


uz 


287, 


,30 


-2. 


.30 


-20. 


,9 


0. 


,7164 


9. 


,331 


0. 


.198 


CAR 


WW 


288. 


.20 


0. 


.00 


-13. 


,0 


0. 


,6700 


9. 


,794 


0. 


.388 


CAR 


wz 


289. 


.20 


-1. 


.20 


-14. 


.7 


1 


,3619 


9. 


.255 


0. 


.376 


CAR 


XX 


291, 


,30 


-4. 


.90 


-10. 


,9 


1, 


,1963 


9. 


.341 


0. 


.344 


CAR 


XY 


291, 


,40 


-3. 


.90 


-5. 


.6 


1. 


,0946 


9. 


.334 


0. 


.405 


CAR 


xz 


290. 


.30 


-0. 


.80 


1. 


,5 


1. 


.2214 


8. 


.595 


0. 


.360 


CAR 


YZ 


285. 


.60 


-1. 


.40 


1. 


,0 


1. 


,2592 


8. 


,709 


0. 


.386 
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Table [| — Continued 
Cepheid £ b v r a log P (V) b > c E(B-V) 



CAR 


AQ 


285. 


,80 


-3 


.30 


2. 


. 1 


0. 


,9899 


8. 


.821 


0. 


, 175 


CAR 


CC 


289. 


.40 


-1 , 


.60 


6. 


. 9 





,6776 


12, 


.021 


0, 


.503 


CAR 


CF 


289. 


,40 


-1 , 


.60 


-3, 


.8 


0. 


,7400 


12, 


.400 


0. 


,602 


CAR 


CN 


283. 


,60 


-1 , 


.30 


9. 


. 


0. 


,6931 


10. 


.670 


0. 


.407 


CAR 


CQ 


286. 


.20 


-1 , 


.70 


20. 


, 6 





.7258 


13. 


.400 


0. 


,902 


CAR 


CR 


285. 


.70 


-0. 


.40 


25. 


,3 


0. 


.9895 


11 . 


.566 


0. 


,484 


CAR 


CS 


288. 


,60 


0. 


.20 


10, 


,8 


0. 


,8236 


12. 


.400 


0. 


.702 


CAR 


CY 


289. 


,50 


-0, 


.90 


8, 


, 6 


0. 


,6300 


9. 


.782 


0. 


,380 


CAR 


DY 


288. 


.80 


-1 . 


.00 


-1 . 


,3 





,6698 


11 . 


.394 


0. 


.382 


CAR 


ER 


290. 


. 10 


1 . 


.50 


-19. 


,0 


0. 


.8876 


6. 


.807 


0. 


. 121 


CAR 


FF 


286. 


.90 


0. 


.60 


-14. 


,5 


1 . 


.2130 


12. 


.300 


0. 


.900 


CAR 


FH 


288. 


,40 


-1 . 


.30 


14, 


, 


0. 


,7531 


12. 


.500 


0. 


.700 


CAR 


FI 


287. 


,80 


0. 


.70 


13. 


, 2 


1 . 


, 1289 


11 . 


.655 


0. 


,683 


CAR 


FK 


288. 


,70 


-0. 


.40 


29, 


, 9 


1 


.3665 


12. 


.200 


0. 


,900 


CAR 


FM 


289. 


,90 


-1 . 


.00 


39, 


,7 


0. 


.8830 


12. 


.400 


0. 


,800 


CAR 


FN 


289. 


,60 


-0. 


. 10 


-13, 


, 


0. 


,6614 


11 . 


.623 


0. 


.553 


CAR 


FO 


290. 


,50 


-2 


. 10 


-10, 


,8 


1 . 


,0152 


10. 


.773 


0. 


,445 


CAR 


FQ 


290. 


.90 


-0. 


.40 


-3. 


. 1 


1 


,0117 


11 . 


.821 


0. 


.811 


CAR 


FR 


291 . 


. 10 


0. 


.60 


-7. 


,3 


1 


,0301 


9. 


.675 


0. 


.346 


CAR 


FZ 


288. 


,40 


0. 


.30 


-1 , 


. 6 


0. 


,5536 


11 . 


.600 


0. 


,619 


CAR 


GI 


290. 


,30 


2. 


.50 


-20, 


.6 


0. 


,6465 


8. 


.320 


0. 


. 188 


CAR 


GS 


289. 


,20 


-1 . 


.80 


32, 


. 1 


0. 


,6080 


12. 


.500 


0. 


,700 


CAR 


GZ 


284. 


.70 


-1 . 


.90 


-8, 


.5 





.6190 


10. 


.239 


0. 


,407 


CAR 


HS 


285. 


.30 


-1 . 


.80 


5. 


, 9 





.7069 


12. 


.300 


0. 


.700 


CAR 


IM 


289. 


, 10 


-0, 


.80 


3. 


.8 


0. 


,7272 


12. 


.300 


0. 


.700 


CAR 


10 


289. 


,30 


-0, 


.90 


39, 


,3 


1 . 


, 1337 


10. 


.938 


0. 


,453 


CAR 


IT 


291. 


,50 


-1. 


.10 


-14. 


.9 





,8770 


8. 


.094 


0. 


.204 


CAS 


RS 


114. 


.50 


0, 


.80 


-24. 


,5 


0. 


.7991 


9. 


.928 


0. 


.818 


CAS 


RW 


129, 


,00 


-4. 


.60 


-71, 


.3 


1. 


,1702 


9. 


.224 


0. 


.408 


CAS 


RY 


115. 


,30 


-3. 


.30 


-70, 


.5 


1. 


,0840 


9. 


.957 


0. 


.614 


CAS 


SU 


133. 


,50 


8. 


.50 


-7, 


,0 


0. 


.2898 


5. 


.970 


0. 


.288 


CAS 


SW 


109, 


,70 


-1. 


.60 


-38, 


,0 


0. 


.7357 


9. 


.700 


0. 


.475 


CAS 


SY 


118. 


,20 


-4. 


.10 


-47, 


.1 


0. 


.6097 


9. 


.858 


0. 


.448 


CAS 


SZ 


134, 


,80 


-1, 


.20 


-45, 


,9 


1. 


,1344 


9. 


.852 


0. 


,767 


CAS 


uz 


125, 


,50 


-1, 


.60 


-51, 


.0 


0. 


,6293 


11. 


.351 


0. 


.498 


CAS 


vv 


130. 


.40 


-2. 


.10 


-50. 


.5 





,7930 


10. 


.751 


0. 


.528 


CAS 


vw 


124. 


.60 


-1, 


.10 


-58. 


.5 





,7777 


10. 


.708 


0. 


.458 


CAS 


XY 


122, 


,80 


-2, 


.80 


-42, 


.0 


0. 


,6534 


9. 


.979 


0. 


.534 


CAS 


AP 


120. 


,90 


0. 


.10 


-44, 


.5 


0. 


,8355 


11. 


,528 


0. 


,781 
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Table [| — Continued 
Cepheid £ b v r a log P (V) b > c E(B-V) 



CAS 


BP 


125. 


,40 


2. 


.80 


-46, 


,5 


0, 


.7974 


10. 


.916 


0. 


,883 


CAS 


CD 


115. 


.50 


1 . 


. 10 


-54, 


,9 


0, 


.8922 


10. 


,746 


0, 


.766 


CAS 


CF 


116. 


,60 


-1 , 


. 10 


-77, 


,7 


0. 


,6880 


11 , 


, 126 


0. 


.539 


CAS 


CG 


116, 


,80 


-1 , 


.30 


-79, 


,3 


0. 


,6401 


11 , 


,355 


0. 


.663 


CAS 


DD 


116, 


.80 


0, 


.50 


-70, 


, 1 


0. 


,9918 


9. 


,878 


0. 


.481 


CAS 


DF 


136, 


.00 


1 , 


.50 


-33. 


,0 


0. 


,5834 


10. 


,856 


0. 


,569 


CAS 


DL 


120, 


,30 


-2 


.60 


-38. 


, 1 


0. 


,9031 


8. 


,964 


0. 


,510 


CAS 


FM 


117. 


,80 


-6, 


.20 


-29. 


, 1 


0. 


,7641 


9. 


, 128 


0. 


,346 


CAS 


V636 


127. 


.50 


1 , 


. 10 


-24. 


,9 


0. 


,9231 


7. 


. 186 


0. 


.786 


CEN 


V 


316. 


.40 


3, 


.30 


-23. 


,2 


0. 


,7399 


6. 


,820 


0. 


.290 


CEN 


TX 


315. 


.20 


-0. 


.60 


-52. 


,0 


1 


,2328 


10. 


.537 


0. 


.982 


CEN 


UZ 


295. 


,00 


-0, 


.90 


-11 . 


,9 


0. 


,5230 


8. 


.760 


0. 


.278 


CEN 


vw 


307, 


,60 


-1 . 


.60 


-30. 


,8 


1 . 


, 1771 


10. 


.242 


0. 


.433 


CEN 


XX 


309, 


,50 


4. 


.60 


-18. 


,8 


1 . 


,0397 


7. 


.799 


0. 


.264 


CEN 


AY 


292. 


,60 


0. 


.40 


-15. 


,5 


0. 


.7251 


8. 


,820 


0. 


.309 


CEN 


AZ 


292. 


,80 


-0, 


.20 


-11 . 


,5 


0. 


,5066 


8. 


,635 


0. 


. 174 


CEN 


BB 


296. 


,40 


-0, 


.70 


-15. 


,8 


0. 


,6017 


10, 


, 146 


0. 


.386 


CEN 


BK 


298. 


.00 


-1 . 


.00 


-26. 


.3 





,5016 


10. 


.063 


0. 


.327 


CEN 


IZ 


294. 


.90 


-0. 


.50 


-19. 


,7 





,7703 


12. 


,500 


0. 


.800 


CEN 


KK 


294. 


,20 


2. 


.70 


-2. 


,7 


1 . 


,0856 


11 , 


.502 


0. 


,608 


CEN 


KN 


307, 


,80 


-2 


. 10 


-39 


,7 


1. 


,5317 


9. 


,855 


0. 


.654 


CEN 


LV 


294, 


,30 


-1 . 


.70 


-16, 


,7 


0. 


,6968 


12, 


. 100 


0. 


.700 


CEN 


MY 


305. 


.30 


1 . 


.20 


-15, 


,2 


0. 


,5704 


12. 


,018 


1 . 


, 158 


CEN 


MZ 


305. 


.40 


-1 . 


.60 


-30, 


, 6 


1 . 


.0151 


11 . 


,546 


0, 


.852 


CEN 


00 


306. 


,90 


-0, 


.60 


-35, 


,7 


1 . 


, 1099 


12, 


.004 


1 . 


. 107 


CEN 


QY 


311. 


,90 


-0, 


.20 


-58, 


,7 


1 . 


,2492 


11 , 


,796 


1 . 


.398 


CEN 


V339 


313. 


.50 


-0, 


.50 


-24, 


.4 





,9762 


8, 


,689 


0. 


.415 


CEN 


V378 


306. 


.10 


0, 


.30 


-16, 


,5 


0. 


,8102 


8, 


,464 


0. 


.386 


CEN 


V381 


310. 


,80 


4. 


.40 


-31. 


,8 


0. 


,7058 


7, 


,659 


0. 


.214 


CEN 


V419 


292, 


,10 


4. 


.30 


-15. 


,2 


0. 


,7410 


8, 


,181 


0. 


.188 


CEP 


DELTA 


105, 


,20 


0. 


.50 


-17. 


,4 


0. 


.7297 


3. 


.954 


0. 


.113 


CEP 


AK 


105, 


,10 


0. 


.40 


-45. 


,6 


0. 


,8593 


11. 


,186 


0. 


,664 


CEP 


CP 


100. 


,40 


1. 


.10 


-40. 


,5 


1. 


,2518 


10. 


,619 


0. 


,644 


CEP 


CR 


107. 


,60 


0. 


.30 


-32. 


,3 


0, 


.7947 


9, 


,629 


0. 


,712 


CIR 


AX 


315. 


,80 


4, 


.00 


-25. 


,5 


0, 


.7221 


5. 


,880 


0. 


,167 


CRU 


R 


299. 


.60 


1. 


.10 


-16. 


,5 


0. 


.7654 


6. 


,771 


0. 


.203 


CRU 


S 


303. 


.30 


4. 


.40 


-7. 


.1 





,6712 


6. 


,599 


0. 


.177 


CRU 


T 


299, 


,40 


0. 


.40 


-9. 


,8 


0. 


,8282 


6. 


,570 


0. 


.204 


CRU 


X 


302, 


,30 


3. 


.80 


-25. 


,0 


0. 


,7938 


8. 


,357 


0. 


.287 
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CRU 


SU 


299. 


,20 


-0. 


,60 


-33, 


, 6 


1 . 


, 1088 


9. 


.822 


0. 


.930 


CRU 


sv 


296. 


.80 


-0. 


.40 


-15. 


.4 





,8453 


12, 


. 130 


0. 


.790 


CRU 


TY 


298. 


,00 


-0. 


.20 


-14, 


,3 


0. 


,6980 


11 . 


.700 


0. 


.800 


CRU 


VV 


299. 


,90 


-1 . 


,80 


-35. 


. 1 


0. 


,7868 


12. 


.300 


1 . 


.301 


CRU 


vw 


300. 


.90 


-0. 


.70 


-2. 


,7 


0. 


.7214 


9. 


.602 


0. 


.638 


CRU 


vx 


300. 


.90 


1 . 


.60 


-28. 


,5 


1 . 


,0868 


12. 


,007 


0. 


,891 


CRU 


AD 


298. 


,50 


0. 


,50 


-34, 


.3 


0. 


.8060 


11 . 


.037 


0. 


,642 


CRU 


AG 


301 , 


,70 


3. 


, 10 


-8, 


, 5 


0. 


,5840 


8. 


.204 


0. 


,220 


CRU 


BG 


300. 


.40 


3. 


.40 


-20. 


,3 


0. 


.5241 


5. 


.462 


0. 


.078 


CYG 


X 


76. 


.90 


-4. 


.30 


8, 


.3 


1 . 


,2145 


6. 


.396 


0. 


.289 


CYG 


SU 


64. 


.80 


2. 


.50 


-21 . 


.4 





,5850 


6. 


.857 


0. 


. 116 


CYG 


SZ 


84. 


,40 


4. 


,00 


-12, 


, 1 


1 . 


, 1792 


9. 


,434 


0. 


.598 


CYG 


TX 


84. 


,40 


-2. 


,30 


-20, 


.4 


1 . 


, 1676 


9. 


.515 


1 . 


.093 


CYG 


VX 


82. 


,20 


-3. 


.50 


-18. 


,0 


1 


.3038 


10. 


.073 


0. 


.742 


CYG 


VY 


82. 


,90 


-4. 


.60 


-11 . 


, 9 


0. 


.8953 


9. 


,590 


0. 


.617 


CYG 


VZ 


91, 


,50 


-8. 


.50 


-18, 


, 5 


0. 


,6871 


8. 


.957 


0. 


.290 


CYG 


BZ 


84. 


,80 


1 . 


.40 


-13. 


, 2 


1 . 


,0061 


10. 


.232 


0. 


.870 


CYG 


CD 


71 . 


. 10 


1 . 


.40 


-11 . 


.6 


1 . 


,2323 


8. 


.957 


0. 


.493 


CYG 


DT 


76, 


.50 


-10 


,80 


-1 . 


. 6 


0. 


,3978 


5. 


.774 


0. 


.065 


CYG 


GH 


66. 


,50 


-0. 


. 10 


-11 , 


. 1 


0. 


,8931 


9. 


.937 


0. 


.626 


CYG 


MW 


70, 


,90 


-0. 


.60 


-16. 


.4 


0. 


,7749 


9. 


.485 


0. 


.642 


CYG 


V386 


85, 


,50 


-4. 


,90 


-5, 


. 1 


0. 


,7207 


9. 


.631 


0. 


.836 


CYG 


V402 


74. 


. 10 


2. 


.30 


-12. 


,7 





,6400 


9. 


.871 


0. 


.405 


CYG 


V438 


77. 


.60 


2. 


.30 


-10. 


, 2 


1 . 


.0496 


10. 


.939 


1 . 


.200 


CYG 


V459 


90. 


,50 


0. 


.70 


-20, 


,7 


0. 


,8604 


10. 


.601 


0. 


.748 


CYG 


V520 


87, 


,50 


1 . 


,60 


-22, 


,8 


0. 


,6073 


10. 


.853 


0. 


.752 


CYG 


V532 


89. 


.00 


-3. 


.00 


-16. 


.2 





,5164 


9. 


.090 


0. 


.511 


CYG 


V1334 


83. 


.60 


-8. 


.00 


-5. 


,2 


0. 


,5226 


5. 


.885 




.037 


CYG 


V1726 


92, 


,51 


-1. 


.61 


-15. 


.3 


0. 


,6269 


9. 


.011 


0. 


.311 


DDR 


BETA 


271. 


,70 


-32. 


,80 


7. 


.3 


0. 


,9931 


3. 


.755 


0. 


.070 


GEM 


ZETA 


195, 


.70 


11. 


.90 


6. 


,8 


1 


.0066 


3. 


.918 


0. 


.046 


GEM 


W 


197. 


,40 


3. 


.40 


-0. 


,1 


0. 


.8985 


6. 


.951 


0. 


.285 


GEM 


RZ 


187. 


.70 


-0. 


.10 


12. 


.0 


0. 


.7427 


10. 


.016 


0. 


.543 


GEM 


AA 


184, 


,60 


2. 


.70 


9. 


.5 


1. 


,0532 


9. 


.735 


0. 


.327 


GEM 


AD 


193, 


,30 


7. 


.60 


45, 


.0 


0. 


,5784 


9. 


.858 


0. 


.180 


GEM 


BB 


199. 


.40 


2. 


.30 


64. 


.3 





,3632 


11. 


.364 


0. 


,444 


GEM 


DX 


198. 


.10 


2. 


.80 


16. 


.0 





,4964 


10. 


.742 


0. 


.436 


LAC 


V 


106, 


,50 


-2. 


.60 


-25. 


.4 


0. 


,6976 


8. 


.939 


0. 


,350 


LAC 


X 


106, 


,50 


-2. 


.50 


-25, 


.3 


0. 


,7359 


8. 


.407 


0. 


.356 
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LAC 


Y 


98 , 


, 70 


-4. 


.00 


-22 , 


, 


. 


, 6359 


9 , 


. 150 


. 


,225 


LAC 


Z 


105 . 


, 80 


-1 . 


.60 


-30. 


.0 


1 


, 0369 


8, 


.418 


. 


.394 


LAC 


RR 


105. 


,60 


-2. 


,00 


-38, 


,8 


0. 


,8073 


8, 


.844 


0. 


,348 


LAC 


BG 


93. 


, 00 


-9 . 


.30 


-18, 


. 6 


0. 


,7269 


8, 


.878 


. 


,332 


LUP 


GH 


324. 


.90 


3. 


.30 


-16. 


, 1 


0. 


.9678 


7. 


,633 


0. 


.358 


MDN 


T 


203 . 


, 60 


-2. 


.60 


22 . 


. 


1 . 


.4317 


6 . 


, 129 


. 


.218 


MDN 


SV 


203. 


,70 


-3. 


.70 


27, 


, 1 


1 . 


, 1828 


8. 


,260 


0. 


,254 


MDN 


TX 


214. 


, 10 


-0. 


.80 


51 , 


, 


. 


, 9396 


10, 


,965 


. 


,490 


MDN 


TY 


213. 


.40 


1 . 


.20 


59. 


. 5 





,6055 


11 . 


,400 


0. 


,700 


MDN 


TZ 


214. 


.00 


1 . 


.30 


34. 


.0 





,8709 


10. 


,775 


0. 


,427 


MDN 


XX 


215 . 


.50 


-1 . 


. 10 


66. 


. 4 





, 7370 


11 . 


.898 


. 


. 566 


MDN 


AA 


217, 


, 10 


-0. 


.30 


66, 


. 1 


0. 


,5953 


13, 


.300 


0. 


.779 


MDN 


BE 


204. 


,90 


1 . 


.00 


45 , 


. 4 


. 


,4322 


10, 


.579 


. 


, 590 


MDN 


CV 


208. 


,60 


-1 . 


.80 


18. 


, 9 


0. 


.7309 


10, 


,306 


0. 


,673 


MDN 


EE 


220 . 


,00 


3 . 


.80 


68. 


, 9 


. 


. 6821 


12. 


,500 


. 


,469 


MDN 


FG 


221 . 


,70 


-0. 


. 10 


88, 


,7 


0. 


,6529 


12, 


,900 


1 . 


,000 


MDN 


FI 


221 . 


,50 


1 . 


.00 


85 , 


, 8 


. 


, 5169 


12 , 


,931 


. 


,515 


MDN 


V465 


214. 


.90 


3. 


.70 


60. 


.7 





,4334 


10. 


,380 


0. 


.260 


MDN 


V508 


208. 


.90 


. 


.90 


58. 


. 8 





, 6164 


10. 


,551 


. 


. 321 


MUS 


R 


302. 


, 10 


-6. 


.50 


0. 


. 1 


0. 


,8756 


6. 


,285 


0. 


, 138 


MUS 


S 


299. 


,60 


-7. 


.50 


-1 , 


.7 


0. 


,9849 


6. 


, 123 


0. 


, 162 


MUS 


RT 


296. 


,50 


-5 . 


.30 


-5 , 


. 5 


. 


,4894 


8. 


,977 


. 


,325 


MUS 


TZ 


296. 


.60 


-3. 


.00 


-24. 


. 6 





.6942 


11 . 


,759 


0. 


.657 


MUS 


UU 


296 . 


. 80 


-3. 


.20 


-17. 


. 


1 . 


. 0658 


9 . 


,772 


. 


.402 


NDR 


S 


327, 


,80 


-5. 


.40 


5. 


.8 


0. 


,9892 


6. 


,430 


0. 


,200 


NDR 


U 


325. 


,60 


-0. 


.20 


-21 , 


. 8 


1 . 


, 1018 


9 . 


,249 


. 


,833 


NDR 


RS 


329. 


.10 


-1. 


.20 


-40. 


.5 





,7923 


10. 


.035 


0. 


.552 


NDR 


SY 


327. 


.50 


-0. 


.70 


-32. 


,9 


1 


,1019 


9. 


,472 


0. 


.745 


NDR 


TW 


330, 


,40 


0. 


.30 


-56, 


.6 


1. 


,0328 


11, 


.670 


1, 


,234 


NDR 


GU 


330, 


,50 


-1, 


.70 


-24, 


,5 


0. 


,5382 


10, 


,425 


0. 


.646 


NDR 


QZ 


329. 


,44 


-02. 


.12 


-38. 


,6 


0. 


.7300 


8, 


,866 


0, 


.278 


NDR 


V340 


329. 


,72 


-02. 


.27 


-40. 


,0 


1. 


.0526 


8. 


,381 


0. 


.332 


OPH 


Y 


20. 


,60 


10. 


.10 


-6. 


,6 


1. 


.2336 


6. 


,167 


0. 


.620 


DPH 


BF 


9. 


,90 


7. 


.10 


-29, 


.0 


0. 


,6094 


7. 


,331 


0. 


.252 


DRI 


RS 


196, 


,60 


0. 


.30 


40, 


,5 


0. 


,8789 


8. 


,412 


0. 


.380 


DRI 


CR 


195, 


.90 


-3. 


.90 


22. 


,4 





,6912 


12. 


,296 


0. 


.538 


DRI 


CS 


198. 


.00 


-4. 


.50 


15. 


.5 


0. 


,5898 


11. 


,390 


0. 


,392 


DRI 


GQ 


199, 


,80 


-4. 


.30 


43, 


,9 


0. 


,9353 


8. 


,965 


0. 


.281 


PER 


SV 


162. 


,60 


-1. 


.50 


-4, 


,5 


1. 


,0465 


8. 


,989 


0. 


,421 
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PER 


SX 


158. 


,90 


-6. 


,40 


5. 


, 5 


0. 


,6325 


11 , 


. 166 


0. 


.471 


PER 


ux 


133. 


.60 


-3. 


. 10 


-41 . 


.5 





,6595 


11 . 


,608 


0. 


.514 


PER 


UY 


135, 


,90 


-1 . 


,40 


-45, 


, 


0. 


,7296 


11 , 


,339 


0. 


.857 


PER 


vx 


132. 


,80 


-3. 


,00 


-35. 


.4 


1 . 


,0372 


9. 


,305 


0. 


.494 


PER 


VY 


135. 


, 10 


-1 . 


.70 


-39. 


,5 


0. 


.7429 


11 . 


,255 


0. 


.877 


PER 


AS 


154. 


. 14 


-0. 


.88 


-25. 


,5 





.6966 


9. 


,726 


0. 


.672 


PER 


AW 


166. 


,60 


-5. 


.40 


6. 


, 9 


0. 


,8104 


7. 


,486 


0. 


.511 


PER 


V440 


135, 


,90 


-5. 


,20 


-26, 


, 1 


0. 


,8791 


6. 


,247 


0. 


.276 


PUP 


X 


236. 


. 10 


-0. 


.80 


65. 


.3 


1 


,4143 


8. 


,490 


0. 


.429 


PUP 


RS 


252. 


.40 


-0. 


.20 


22. 


, 1 


1 


,6169 


7. 


,028 


0. 


.431 


PUP 


VW 


235. 


.36 


-00 


,62 


+24, 


,0 


0. 


,6320 


11 . 


,382 


0. 


.493 


PUP 


VX 


237, 


,00 


-1 . 


.30 


8. 


.8 


0. 


,4789 


8, 


.315 


0. 


. 152 


PUP 


VZ 


243. 


,40 


-3. 


.30 


63, 


,3 


1 . 


,3648 


9. 


.595 


0. 


.454 


PUP 


wx 


241 . 


,20 


-1 . 


.40 


54, 


, 6 


0. 


.9512 


9. 


,058 


0. 


.317 


PUP 


wz 


241 . 


,80 


3, 


.30 


64, 


,0 





.7013 


10. 


,301 


0. 


,228 


PUP 


AQ 


246. 


,20 


0. 


.00 


59, 


,5 


1 . 


,4774 


8. 


.755 


0. 


.491 


PUP 


AT 


254. 


,30 


-1 . 


,60 


26, 


,7 


0. 


,8238 


8. 


.001 


0. 


. 195 


PUP 


BM 


244. 


.50 


-1 . 


.00 


64. 


.4 





,8573 


10. 


.811 


0. 


.581 


PUP 


BN 


247. 


.90 


1 , 


.00 


65. 


,7 


1 . 


, 1359 


9. 


.922 


0. 


.424 


PUP 


HW 


244. 


,80 


0. 


.80 


116, 


.2 


1 . 


, 1288 


12, 


. 129 


0. 


,681 


PUP 


LS 


246. 


,40 


0. 


. 10 


77, 


.4 


1 . 


, 1506 


10, 


.363 


0. 


,460 


PUP 


MY 


261 . 


,30 


-12. 


,90 


12, 


.7 


0. 


,7555 


5. 


.666 


0. 


,088 


SGE 


S 


55. 


.20 


-6. 


. 10 


-10. 


, 1 


0. 


.9233 


5. 


,622 


0. 


. 144 


SGE 


GY 


54. 


.90 


-0. 


.60 


15, 


, 6 


1 . 


,7107 


10, 


,208 


1 . 


,268 


SGR 


U 


13. 


,70 


-4. 


,50 


2. 


, 6 


0. 


,8290 


6. 


,714 


0. 


.393 


SGR 


W 


1 . 


.60 


-4. 


.00 


-28, 


,3 


0. 


,8805 


4. 


,669 


0. 


. 130 


SGR 


X 


1. 


,20 


0. 


.20 


-13. 


,7 


0. 


.8458 


4. 


,561 


0. 


.207 


SGR 


Y 


12. 


.80 


-2. 


.10 


-2. 


,5 


0. 


.7614 


5. 


,742 


0. 


.214 


SGR 


VY 


10. 


,13 


-01. 


,07 


-6, 


,0 


1. 


,1322 


11, 


,529 


1. 


.185 


SGR 


WZ 


12. 


,10 


-1. 


,30 


-15, 


,7 


1. 


,3395 


8, 


,030 


0. 


.450 


SGR 


XX 


15. 


,00 


-1. 


.90 


2. 


.0 


0. 


.8078 


8, 


,851 


0. 


.519 


SGR 


YZ 


17. 


,80 


-7. 


.10 


18. 


.5 


0. 


.9742 


7. 


,350 


0. 


.293 


SGR 


AP 


8, 


,10 


-2. 


.40 


-15. 


,0 


0. 


.7040 


6. 


.942 


0. 


.203 


SGR 


AV 


7. 


,50 


-0. 


.60 


0. 


.0 


1. 


,1878 


11, 


.425 


1. 


.170 


SGR 


AY 


13. 


,30 


-2. 


.40 


-26, 


,5 


0. 


,8176 


10, 


.520 


0. 


.857 


SGR 


BB 


14. 


.70 


-9. 


.00 


7. 


.6 
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a See text for references 

b Caldwell & Coulson (1987) and references therein 
c Pont et al. (1994) 
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Fig. 1. — The pixel location of the center of the spectrum, along the cross-dispersed direction, 
plotted as a function of telescope hour angle. A significant shift is seen near the meridian (with 
some hysteresis), possibly caused by motion of the grism. 
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Fig. 2. — Velocities of the template calculated from spectra of radial velocity standards. The open 
and filled points correspond to velocities measured from the two different locations of the spectrum 
on the chip; the error of each measurement is typically 2.7 km s -1 . The solid line shows the adopted 
mean velocity for the template. 
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Fig. 3. — Radial velocity curves for the newly discovered Cepheids, using the designations of 
Caldwell et al. The curve shape was determined from the period, and its position in phase and 7 
velocity were fit to the radial velocities shown. 
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Fig. 4. — A comparison of the reddenings from Fernie (1990) and (a) Caldwell & Coulson (1987); 
(b) Laney & Stobie 1994. A linear fit to each is shown, revealing a reddening scale difference in 
both. 
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Fig. 5. — A schematic of Milky Way rotation, with labels indicating quantities discussed in the text. 
The sun is indicated with a circle near the top of the figure, and a fiducial Cepheid is indicated 
with a star. The Galactic center is labeled GC. 
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Fig. 6. — Locations of Milky Way Cepheids used in the rotation curve models. Squares indicate 
Cepheids modeled by CC. Open circles indicate additional Cepheids with new reddenings, filled 
circles Cepheids with new radial velocities from MCS and PMB, and stars indicate Cepheids newly 
discovered by Caldwell et al. Cartesian coordinates are shown in units of Rq with the Galactic 
center at (0,0) and the Sun at (0,1). The solar circle (r = 1) is shown with a dotted line. 
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Fig. 7. — Constant b\ 2 contours for two models projected in the 2AR$, mo plane. Contours shown 
are (from inside to outside) l-cr, 90%, 2-cr, 99%, 99.9%, and 99.99%. The top panel is a model using 
all Cepheids with F-band data. The bottom panel is a model that also includes CKS Cepheids, 
where distances were determined from K-band data. 
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Fig. 8. — A plot of the best-fit model parameters as a function of the symmetric ellipticity 
component c(Rq). The left panel shows the fit distance modulus of the galactic center, (m — M)q, 
the right panel circular velocity. A flat rotation curve and constant ellipticity are assumed. The two 
horizontal lines represent 1-a errors on other, non- kinematic distance measurements (Reid 1993). 
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Fig. 9. — Residual velocities of the stars in Model Bl, smoothed with a 0.5 kpc gaussian filter. 
Contour levels are at 2 km s _1 increments. Shaded contours are negative velocity residuals, i.e. the 
mean motion of stars in these areas is toward the sun with respect to the model. 
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Fig. 10. — A schematic of the closed elliptical orbits that pass through the sun in our fiducial 
model. The general orientation of the orbits is depicted; the ellipticity has been exaggerated for 
clarity. 



